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Abstract 

In this work we report on the effects of short-range correlations 
upon the matrix elements of neutrinoless double beta decay (Oi^/3/3). 
We focus on the calculation of the matrix elements of the neutrino- 
mass mode of Qv(3f3 decays of ^^Ca and ^^Ge. The nuclear-structure 
components of the calculation, that is the participant nuclear wave 
functions, have been calculated in the shell-model scheme for ^'^Ca 
and in the proton-neutron quasiparticle random-phase approximation 
(pnQRPA) scheme for ^^Ge. We compare the traditional approach of 
using the Jastrow correlation function with the more complete scheme 
of the unitary correlation operator method (UCOM). Our results in- 
dicate that the Jastrow method vastly exaggerates the effects of short- 
range correlations on the Qv(3(3 nuclear matrix elements. 



PACS number(s): 21.60.Cs, 23.40.Hc, 27.40.+Z. 

Keywords: Neutrinoless double beta decay, nuclear matrix elements, short- 
range correlations, unitary correlation operator method 

To appear in Physics Letters B (2007) 



A good knowledge of the nuclear matrix elements governing the decay 
rates of neutrinoless double beta (Oz//3/3) decay is mandatory if one wants to 
extract information about the neutrino mass from the current experimental 
limits for the half-lives of 0z//9/9-decay transitions [1, 2]. The standard theo- 
retical methods which are suitable for the calculation of the relevant nuclear 
matrix elements can be found in the literature, e.g. in [3, 4]. Although many 
of the formalisms are well established, difficulties arise from the approxima- 
tions which are needed in order to perform the actual calculations. 

In the mass mode of the Ou^P decay the involved two nucleons exchange a 
light Majorana neutrino [5]. Average value of the exchanged momentum is of 
the order of 100 — 200 McV/c and thus the involved nucleons arc on average 
at close distance from each other. There is, however, a minimum relative 
distance of the order of 1 fm after which the two nucleons may eventually 
overlap. In nuclear matter this overlapping cannot happen and in theoretical 
description of the 0z//9/9 decay one needs to take into account this fact. Based 
on this it has been argued [6] that special measures have to be taken when 
performing nuclear-structure calculations using the mean-field picture with 
residual two-body interactions between the two interacting nucleons. In the 
case of the OuPP decay these measures boil down to introducing an explicit 
Jastrow type of correlation function into the involved two-body transition 
matrix elements [3]. Using this method in the numerical calculations of 
Oi//3/3-decay matrix elements considerable corrections to the involved Fermi 
and Gamow-Teller nuclear matrix elements were reported [3, 7]. 

In [8] a different method was used to explicitly take into account the short- 
range correlations. This approach is based on the use of the Horie-Sasaki 
method to evaluate the involved radial form factors and the short-range cor- 
relations were considered to arise from the oo exchange in the nucleon-nucleon 
interaction [9, 10]. Contrary to [3, 7], in this approach only relatively small 
corrections to the involved nuclear matrix elements were obtained [8, 11]. 
Instead of using the above described methods, one can use the more com- 
plete new concept of unitary correlation operator method (UCOM) [12] to 
take into account the short-range effects in Ou^P decay. In this method a 
unitary correlation operator moves a pair of nucleons away from each other 
whenever they start to overlap. This method also conserves the probability 
normalization of the relative wave function. 

In this work we address the important issue of short-range correlations 
in the computation of nuclear matrix elements involved in the neutrinoless 
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double beta decay. We have used both the UCOM and Jastrow methods and 
we compare the results for the ground-state-to-ground-state OuPP decays of 
^^Ca and ^^Ge. For ^^Ca we calculate the relevant nuclear wave functions in 
the solid theoretical framework of the nuclear shell model. In order to accom- 
plish this we have used the OXBASH code [14], which is actually available 
to any practitioner in the field. For ^^Gc we have used the framework of 
the proton-neutron quasiparticle random-phase approximation (pnQRPA), 
suitable for calculations of nuclear properties of medium-heavy and heavy 
nuclei. 

Our calculated results show that the reduction caused by the inclusion 
of the short range correlations depends on the multipole which contributes 
to the Oi^PP matrix element. In addition, the strength distributions of the 
multipoles are practically unaltered by the short-range correlations, suggest- 
ing the effect to be just an overall multipole-dependent scaling. This scaling 
factor is close to unity for all multipoles in the UCOM scheme, but acquires 
strongly reduced values for high multipoles in case of the Jastrow method. 

We start the quantitative scrutiny of the effects of short-range correlations 
by briefly presenting the central ideas behind our computations. By assuming 
the neutrino mass mechanism to be the dominant one in the Ou^P decay one 
can write as a good approximation the inverse of the half-life as [4] 

2 / 2 \ ^ 

= G^l^i^j , (1) 

where rrie is the mass of the electron and 

(m.) = EAfm,|C/e,f (2) 
j 

is the effective mass of the neutrino, A^^ being the GP phase. Furthermore, 

the quantity Cf'^^ of Eq. (1) is the leptonic phase-space factor defined in [4]. 
The double Gamow-Teller and double Fermi nuclear matrix-elements in (1) 
are defined as [5] 

Mt^ = j:mh4rmn,Eam), (3) 

a 

= E(0;il^+(^mn,£^a)K-^n)||0+). (4) 
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Here the summation runs over all states a of the intermediate nucleus, which 
in this case are ^^Sc and ''^As. The definition of the neutrino potential 
h+{rmn: Eg) Can be found in Refs. [3, 4, 5]. 

The traditional way [3] to include the short-range correlations in the Oi'(3(3 
matrix elements is by introducing the Jastrow correlation function f]{r). 
It has to be noted that this particular variant of the Jastrow function is a 
rudimentary one and does not do full justice to the name Jastrow correlations. 
For example, in fight nuclei accurate Monte-Carlo calculations are based on 
Jatrow-like correlations wfiich arc variationally determined and have different 
ansatz functions in different isospin channels. This fact notwithstanding we 
choose to call fiere tfie rudimentary approacfi of [3] as Jastrow metfiod since 
tfiis is tfie term adopted by tfie double-beta-decay community. 

Tfie Jastrow function depends on tfie relative distance r = |ri — of two 
nucleons, and in tfie Jastrow scheme one replaces tfie bare Ov^P operator O 
by a correlated operator Oj by tfie simple procedure 

(0+||O||0+) ^ (O;||0a||O+) = (0+||/aO/j||0+). (5) 

A typical cfioice for tfie function /j is 

Mr) = 1 - e— (l - br') , (6) 

witfi a = 1.1 fm^ and b = 0.68 fm^. As a result, tfie Jastrow function effec- 
tively cuts out tfie small r part from tfie relative wave function of tfie two 
nucleons. For tfiis reason, tfie traditionally adopted Jastrow procedure does 
not conserve tfie norm of tfie relative wave function and one sfiould use, in 
principle, tfie operator 

wfien tfie initial and final states are tfie same, fiere denoted by For 
different initial and final states normalization is even more problematic. Even 
a proper normalization does not guarantee tfiat tfie sfiort-range correlations 
would be correctly treated by tfie Jastrow procedure. 

To circumvent tfie difficulties associated witfi tfie use of a Jastrow function 
one can adopt tfie more refined unitary correlation operator metfiod (UCOM) 
[12] wfien discussing sfiort-range effects in tfie OuP(3 decay. In tfie UCOM 
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one obtains the correlated many-particle state |^) from the uncorrelated one 

as 
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where C is the unitary correlation operator. The operator C is a product 
of two unitary operators: C = C^Cr, where describes short-range tensor 
correlations and Cr central correlations. Due to unitarity of the operator C, 
the norm of the correlated state is conserved. Moreover, one finds for the 
matrix element of an operator A 



so that it is equivalent to use either correlated states or correlated operators. 

The exact form of the operator C is gained by finding the minimum of 
the Hamiltonian matrix element (^|C^i/C|\l'). Therefore, the choice of the 
two-body interaction in H affects also the form of C. Explicit expressions 
for the operators Cr and Cn can be found in Refs. [12, 13]. Application 
of these expressions to the double Gamow- Teller and Fermi matrix elements 
shows that the tensor correlations of Cq vanish and we are left with only the 
central correlations. 

For ''^Ca the nuclear-structure calculations were handled by the shell- 
model code OXBASH [14]. In our calculations we have used the FPBP 
two-body interaction of [15], which was obtained by fitting the Kuo-Brown 
interaction to experimental data. Due to the fact that we have limited our 
model space to the pf shell, the Oi'f3f3 matrix elements are composed of only 
positive-parity states. The shell-model calculations had to be truncated by 
requiring that the minimum number of particles in the 0f7/2 orbital be 4. 

Our main results for ^^Ca are presented in Tables 1 and 2. In these tables 
we list the calculated multipole decomposition and total values of the matrix 
elements Mq^^ and Mp'^^ for four different cases. In the first case, which we 
refer to as bare matrix elements, we have not taken into account any short- 
range correlations. In the second case the short-range effects were handled by 
the use of the Jastrow function (6) and the replacement (5). In the third and 
fourth cases we have used the UCOM to account for the short-range effects. 
The Kuo-Brown interaction was not derived via UCOM, as it should be if it 
were to be used in the same calculation as the UCOM-derived double-beta 
operator. To access the magnitude of the resulting effect, we have adopted 
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Table 1: Multipole decomposition and the total value of the matrix clement 
M^°^ for 48 Ca. The cases are: no short-range correlations included (bare), 
with Jastrow correlations and with UCOM correlations using Bonn-A and 
Argonne V18 parametrizations. 

UCOM 



J" 


bare 


Jastrow 


Bonn-A 


AVIS 


1+ 


-0.330 


-0.305 


-0.322 


-0.319 


2+ 


-0.117 


-0.092 


-0.108 


-0.104 


3+ 


-0.327 


-0.246 


-0.302 


-0.293 


4+ 


-0.066 


-0.035 


-0.054 


-0.051 


5+ 


-0.246 


-0.121 


-0.212 


-0.199 


6+ 


-0.042 


-0.008 


-0.030 


-0.027 


7+ 


-0.150 


-0.029 


-0.120 


-0.107 


sum 


-1.278 


-0.835 


-1.150 


-1.101 



two different UCOM parameter sets in the present calculation. These two 
parameter sets were obtained by minimizing the energy for the Bonn-A and 
Argonne VIS potentials. Both of the used UCOM parameter sets can be 
found in [16]. 

As the results in Tables 1 and 2 indicate, the differences between the 
results obtained by the use of the two UCOM parameter sets are small. 
Therefore, we expect that the results obtained by the use of the Kuo-Brown 
UCOM parameters do not deviate significantly from the Bonn-A or Argonne 
V18 results. We also note that there exist a small UCOM contribution to 
the double Fermi matrix element Mp'^^ coming from the odd- J intermediate 
states. This is explained by the fact that in Ref. [16] slightly different 
parameters were given to the S* = and S* = 1 channels. 

In Fig. 1 we show graphically the multipole decomposition of the total 
matrix element Mq^^ — [g^^ / g^f M^^^ of (1) for the four different cases of 
Tables 1 and 2. The ratio gx/gY — —1.254 was used in this plot. As can be 
seen, the results obtained by using the two different UCOM parameter sets 
do not differ significantly. Also, one can see that the effects of the Jastrow or 
UCOM correlations grow with increasing J of the intermediate states. For 
the extreme case of the 1'^ contributions the switching on of the Jastrow 



5 



Table 2: The same as Table 1 but for M^' 



UCOM 





bare 


Jastrow 


Bonn-A 


AV18 


1+ 


0.000 


0.000 


0.000 


0.000 


2+ 


0.185 


0.145 


0.174 


0.169 


3+ 


0.000 


0.000 


-0.001 


-0.001 


4+ 


0.116 


0.061 


0.102 


0.096 


5+ 


0.000 


0.000 


-0.002 


-0.002 


6+ 


0.061 


0.012 


0.050 


0.045 


7+ 


0.000 


0.000 


-0.002 


-0.002 


sum 


0.367 


0.221 


0.324 


0.308 



correlations changes the value of the matrix element Mq^ (7"'") from —0.150 
to —0.029, roughly corresponding to a factor of 5 reduction. At the same time 
the UCOM correlations produce only a 20% - 30% reduction from the bare 
matrix element. It seems that in a situation like this blind use of Jastrow 
correlations cuts out relevant parts of the nuclear many-body wave function. 
From the tables one deduces that the Jastrow correlations cause some 35% - 
40% reduction to the magnitudes of the total matrix elements, whereas the 
UCOM causes a reduction of 10% - 16%. It is worth pointing out that our 
numbers for the Jastrow case coincide with the numbers of the corresponding 
earlier calculation performed by the Strasbourg group [17]. 

To trace the source of differences between the Jastrow and UCOM cor- 
rected matrix elements we show for the "^^Ca decay in Fig. 2 the radial depen- 
dence of the two-particle Gamow-Teller OuPP matrix element in the special 
case oi p — p' — n — n' — O/7/2 and J — 7 (this is the contribution to 
Eq. (4.16) of [3] without including the one-body transition densities and the 
overlap of the two complete sets of pnQRPA states). The oscillator param- 
eter value b = 2.0 fm was used in the plot. For the case of UCOM contri- 
bution we have used the correlated wave functions and the approximation 
r — R-{r) R+ir) — r for illustrative purpose [18]. Thus, the UCOM plot 
should be taken only as a schematic one. From the figure one can see that 
the Jastrow correlations cut out a significant part of the matrix element at 
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Figure 1: Multipole decomposition of the total Oi//3/3 decay matrix element 
- {gy/gpf M^°") for ^^Ca. The cases are: no short-range correlations 
included (bare), with Jastrow correlations and with UCOM correlations using 
the Bonn-A and Argonne V18 parametrizations. 

small r. This leads to a situation, where the total integrated areas under 
the radial curve almost cancel out. In the case of UCOM correlations the 
cancellation is not as severe due to the fact that not so much amplitude is 
lost for small r. 

Our results for the Ot//3/3 decay of ^^Ge are summarized in Fig. 3 and Ta- 
ble 3. The results have been obtained by using the framework of the proton- 
neutron quasiparticle random-phase approximation (pnQRPA) [4, 19]. The 
related calculations, including the BCS and the pnQRPA calculations for 
''^Ge and ™Se, were done in the model space consisting of the single-particle 
lp-0f-2s-ld-0g-0hii/2 orbitals, both for protons and neutrons. The single- 
particle energies were obtained from a spherical Woods-Saxon potential. 
Slight adjustment was done for some of the energies at the vicinity of the 
proton and neutron Fermi surfaces to reproduce better the low-energy spec- 
tra of the neighboring odd-A nuclei and the low-energy spectrum of ''^As. 
The Bonn-A G-matrix [20] was used as a two-body interaction and it was 
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Figure 2: Radial dependence of the two-particle Gamow- Teller Oz//?/? matrix 
element for p = p' = n = n' = O/7/2 and J = 7 in the case of ^^Ca de- 
cay. Shown are the bare matrix element and Jastrow and UCOM correlated 
matrix elements. 

renormalized in the standard way, as discussed e.g. in Refs. [21, 22, 23]. Due 
to this phenomenological renormalization we did not perform an additional 
UCOM renormalization of the two-body interaction. In the present calcu- 
lations we have used the 'default value' g^p = 1.0 for the particle-particle 
interaction parameter of the pnQRPA. 

In Fig. 3 we display for ^^Ge decay the multipole decomposition of the 
total 0z//3/3 matrix element Mq^^ — (gy/ qa)'^ Mp"^ as derived from the pn- 
QRPA calculations. The used symbols are the same as in Fig. 1. The ratio 
Qv/qa = —1.254 was used in the calculations. Since the nuclear wave func- 
tions have been calculated by the use of the Bonn potential, we have used 
only the Bonn-A parametrization for the UCOM. Here one can see a pat- 
tern similar to the case of '^^Ca: the effect of the Jastrow correlations grows 
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Figure 3: The same as Fig. 1 for the decay of ^^Ge calculated by using the 
pnQRPA. Only the Bonn-A parametrization has been used for UCOM. 



strongly with increasing value of the angular momentum of the intermedi- 
ate states. As in the case of the ^^Ca decay the effect is the largest for the 
unnatural-parity states l"*", 2~, S"*", 4~, . . . in an odd-odd nucleus. Contrary 
to the Jastrow-corrected multipole contributions, the UCOM-corrected ones 
stay close to the bare contributions for all intermediate multipoles J^. 

We summarize our results on the 01/(3(3 matrix elements of the ^^Ge decay 
in Table 3, where we give the bare, Jastrow-corrected and UGOM-corrected 
Gamow-Teller, Fermi and total matrix elements. For the total matrix ele- 
ment the Jastrow corrections amount to 30% reduction from the bare matrix 
element, whereas the UCOM corrections are some 7%. This, again, suggests 
that in the earlier calculations [3, 7] the effect of the short-range correlations 
has been considerably overestimated. 

In this Letter we have addressed the important issue of short-range corre- 
lations in the context of neutrinoless double beta decay. We have calculated 
the related nuclear matrix by the nuclear shell model for "'^Ca and by the 
pnQRPA for ^^Ge. The short-range correlations have been calculated by the 
use of the simple Jastrow function and the more refined UCOM method. Our 
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Table 3: Gamow- Teller (M^^^), Fermi (M^^^) and total matrix elements de- 
rived from pnQRPA calculations for the Oi'f3f3 decay of ^^Ge. The cases 'bare', 
'Jastrow' and 'UCOM' are as in Table 1. Only the Bonn-A parametrization 
has been used for UCOM. 

bare Jastrow UCOM Bonn-A 
M^T^ -6.755 -4.681 -6.265 
M^^'"'^ 2.474 1.778 2.310 
total -8.328 -5.811 -7.734 



computed results indicate that the Jastrow method cuts off relevant parts of 
the many-body wave function for high values of angular momentum of the 
intermediate states. This leads to the excessive reduction of 30% - 40% in 
the magnitudes of the nuclear matrix elements. At the same time the UCOM 
reduces the magnitudes of the matrix elements only by 7% - 16%, roughly 
equally for all multipoles. Our results put to question the recent calculations 
where short-range and tensor correlations cause large effects on the nuclear 
matrix elements of neutrinoless double beta decay [7]. Study of the effects 
of the UCOM procedure upon heavier nuclei is in progress. 

Acknowledgements: This work has been partially supported by the 
National Research Council (CONICET) of Argentina and by the Academy 
of Finland under the Finnish Centre of Excellence Programme 2006-2011 
(Nuclear and Accelerator Based Programme at JYFL). We thank also the 
EU ILIAS project under the contract R113-CT-2004-506222. One of the 
authors (O.C.) gratefully thanks for the warm hospitality extended to him 
at the Department of Physics of the University of Jyvaskyla, Finland. 

References 

[1] V. Tretyak, Y. Zdesenko, At. Data Nucl. Data Tables 80 (2002) 83. 
[2] S.R. EUiott, P. Vogel, Ann. Rev. Nucl. Part. Sci. 52 (2002) 115. 
[3] T. Tomoda, Rep. Prog. Phys. 54 (1991) 53. 



10 



[4] J. Suhonen, O. Civitarese, Phys. Rep. 300 (1998) 123. 

[5] M. Doi, T. Kotani, E. Takasugi, Prog. Theor. Phys. Suppl. 83 (1985) 1. 

[6] G.A. Miller, J.E. Spencer, Ann. Phys. 100 (1976) 562. 

[7] V.A. Rodin, A. Faessler, F. Simkovic, P. Vogel, Nucl. Phys. A 766 (2006) 
107. 

[8] F. Krmpotic, J. Hirsch, H. Dias, Nucl. Phys. A 542 (1992) 85. 

[9] G.E. Brown, S.O. Backman, E. Oset, W. Weise, Nucl. Phys. A 286 
(1977) 191. 

[10] IS. Towner, Phys. Rep. 155 (1987) 263. 

[11] C. Barbero, F. Krmpotic, A. Mariano, D. Tadic, Nucl. Phys. A 650 
(1999) 485. 

[12] H. Feldmeier, T. Neff, R. Roth, J. Schnack, Nucl. Phys. A 632 (1998) 
61. 

[13] R. Roth, H. Hergert, P. Papakonstantinou, T. Neff, H. Feldmeier, Phys. 
Rev. C 72 (2005) 034002. 

[14] B.A. Brown, A. Etchegoyen, W.D.M. Rae, The computer code 
OXBASH, MSU-NSCL Report 524 (1988). 

[15] W.A. Richter, M.G. Van Der Merwe, R.E. Juhes, B.A. Brown, Nucl. 
Phys. A 523 (1991) 325. 

[16] T. Neff, H. Feldmeier, Nucl. Phys. A 713 (2003) 311. 

[17] J. Retamosa, E. Caurier, F. Nowacki, Phys. Rev. C 51 (1995) 371. 

[18] In all other numerical applications of the UCOM we have used correlated 
operators without involving any approximations. 

[19] J. Suhonen, From Nucleons to Nucleus: Concepts of Microscopic Nuclear 
Theory, Springer, Berlin, 2007. 

[20] K. Holinde, Phys. Rep. 68 (1981) 121. 

11 



[21] J. Suhonen, T. Taigel, A. Faessler, Nucl. Phys. A 486 (1988) 91. 
[22] J. Suhonen, Nucl. Phys. A 563 (1993) 205. 
[23] J. Suhonen, Nucl. Phys. A 700 (2002) 649. 



12 



